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An analytical study of the behavior of the voter model on the small- world topology is performed. 
In order to solve the equations for the dynamics, we consider an annealed version of the Watts- 
Strogatz (WS) network, where long-range connections are randomly chosen at each time step. The 
resulting dynamics is as rich as on the original WS network. A temporal scale r separates a quasi- 
stationary disordered state with coexisting domains from a fully ordered frozen configuration, r 
is proportional to the number of nodes in the network, so that the system remains asymptotically 
disordered in the thermodynamic limit. 

PACS numbers: 87.23.Ge,89.75.-k,05.70.Ln 



I. INTRODUCTION 



The relationship between nontrivial interaction topologies and ordering phenomena is still a largely unexplored 
topic. The recent burst of activity on complex networks has revealed that many technological, social and biological 
systems have interaction patterns markedly different from structures traditionally studied as regular lattices and 
random graphs 0, 0| • The interest is now focusing on how such complex topologies affect dynamical processes taking 
place on them. Models of ordering dynamics play an important role in this context, since they are commonly used to 
study social phenomena, like cultural assimilation and opinion dynamics |2n^ <S] ; for which the interaction patterns 
are more plausibly described by complex networks than by regular lattices. 

Some ordering processes on complex networks have recently been considered, including the zero temperature Glauber 
dynamics of the Ising model on the Watts-Strogatz network and the Axelrod model on small-world and scale-free 
networks 

In a recent paper |^ , we have studied numerically the dynamics of the voter model on a small- world network. This 
structure, more precisely the Watts-Strogatz (WS) network 0, is one of the simplest examples of complex topology. 
Depending on the parameter p (to be specified below) it interpolates between a one-dimensional lattice with periodic 
boundary conditions (for p = 0) and a random graph (for p = 1). It has been shown |10|. that in a well defined 
range of intermediate values of p, the network has, simultaneously, global properties typical of random graphs (small 
average distance between nodes) and local properties (clustering) typical of regular structures. 

The voter model is possibly the simplest model of an ordering process each site a discrete variable a is 

defined, that may assume two values (tr — ±1) representing two opposite options, for instance the electoral choice in 
favor of two different candidates. Starting from a disordered initial condition, the model follows a simple dynamical 
evolution: at each time step one site is selected at random and set equal to one of its nearest neighbors, chosen at 
random in its turn. On regular lattices, in d = 1 and d — 2 the model converges to an ordered state with all variables 
having the same value, whereas for d > 3 the system reaches a disordered stationary state The voter model 

on complete graphs has been considered recently |3 ■ 

In Ref. we have found that the nontrivial connectivity pattern of the WS network has a deep impact on the 
ordering dynamics of the voter model. In particular, after an initial transient, the system settles in a quasi-stationary 
state with coexisting domains. If the system size is infinite this state persists forever: at odds with naive expectations, 
long-range connections prevent complete ordering from being reached. If the system size is kept finite instead, the 
stationary state has a finite lifetime and the fully ordered state is quickly reached at the end of it. Interestingly, the 
dependence of the lifetime on the system size is such that the ordered state is reached earlier than on a one-dimensional 
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lattice of the same size. This partially restores the intuitive picture that long-range connections should speed up the 
ordering process. 

In this work we analyze the same problem from the analytical point of view. The very simple form of the transition 
rates for the voter model results in equations of motion for the correlation functions that are not coupled with each 
other. The only difficulty is to carry out the average over trajectories for a fixed realization of the network and 
to average over the topology afterward. To overcome this problem we use an annealed approximation, consisting 
in averaging over the topology before averaging over the trajectories. This approximation is exact for an effective 
network, which we call annealed WS network, where long-range interactions are not quenched from the beginning but 
are extracted randomly at each time step. The dynamics of the voter model on the annealed WS network is then 
solved, revealing a phenomenology as rich as in the quenched case. By comparing the exact results on the annealed 
WS network with the numerical simulations on the quenched one, it turns out that the only discrepancy between the 
two cases is the dependence on the parameter p of the correlation length in the stationary disordered state. 

The paper is organized as follows. In the next section we define precisely the voter model dynamics and the WS 
topology on which the process occurs. Section ITTTl is devoted to the analytical solution of the dynamics. The equation 
of motion for the correlation function is derived, the annealed approximation is introduced and the behavior of the 
system is studied. In section Hvl we check numerically the analytical results for the annealed case and compare them 
with simulations of the quenched case. In section we extend the results to the case where each site is initially 
connected to neighbors. The final section contains a short discussion of the findings. 



II. THE MODEL 



We consider a small- world network defined as the superposition of a one-dimensional lattice of L sites with periodic 
boundary conditions and a random graph In general one can start from a lattice with each site linked to v 

neighbors on the right and v on the left. We now consider v — 1, deferring the discussion of generic values of v 
to Section More precisely, site i is initially connected with sites i — \ and i + 1. Then a link is added between 
any pair of non-nearest neighbor sites with probability p/L. In this way the total number of edges in the system 
is i + L{L — 3)/2 • p/L, so that the average degree per site is finite (2 + p) in the thermodynamic limit L oo. 
The generalization to an initial lattice with connections to k nearest neighbors is straightforward. This topology 
slightly differs from the one originally introduced by Watts and Strogatz ,10.] but has the same properties and is more 
amenable to analytical treatment. The topology is fully specified by the adjacency matrix Q{i,j), which is 1 if i and 
j are connected and otherwise. The probability distribution of its elements is 

{SQ{i,j),o for i = j 

5Q{^,,)^ for*=j±l (1) 

f '^Q(»,i),i + (1 - L)^Q(t,j)fl otherwise. 

The voter model dynamics is defined by the transition rates. If we call {u} the spin configuration of the system, 
that is {(t} = {cti, (72, . . . , Ci, . • . , cl} and we indicate with {cr'}i the same configuration with the i-th spin flipped, 
the transition rate from state {cr} to state {cr'li is, in complete analogy with the definition for regular lattices, 

^{{o} - Wh) = ^ (^1 - 1 E ' (2) 

where Zi, the degree of site i, is 

^3 ^^Q{i^j)- 



III. ANALYTICAL TREATMENT 
A. Equation for the correlation function 

Given the explicit expression jSJ of the transition rates, from the master equation 

|p({a}, t) = Y. ^({^'}'^ ^ {^})nwh^ t)~J2 ^({^} ^ Wh)PiW}, t), (3) 

i i 
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one can derive |171 Il8| the equation of motion for the mean spin at site j, Sj = {<Tj 



dt 



and for the two-point correlation function ^ — {(JjCTk) 



E 



k)sk 



dCj,k „^ , Q{i,i)Ck., 
^--2Cm- + ^ — - 



E 



Q{k,i)Cj, 

Zk 



(4) 



(5) 



For p — 0, Eqs. I@J and Q coincide with the equations for the one-dimensional voter model [l^. 

In order to study the dynamics of the voter model on the Watts-Strogatz network we must average Eqs. Q and ^ 
over the disordered topology. Indicating with an overbar averaged quantities A = J Yii j (^Qihj)^ ' P[Qihj)]^ the 
equation for the mean spin is 



dt ^ ^ 



/ k)sk 



and the one for the pair correlation function is 



dCj^k 
dt 



E 



V Zk 



(6) 



(7) 



To evaluate the average values appearing on the r. h. s. of Eqs. © and (O we introduce the annealed approximation 

(8) 



A = 



Q{j,k)\ ^ 



This approximation can be seen as considering annealed transition rates 

1 



i-E 



fQ{i,k) 



V Zi 



(9) 



The evaluation of (^^^^^ is readily performed 



l.rn 



dQ{3,k)P[Q{],k)]Q{],k) [ [] dQ{l,m)P[Q{l,m)] ^ 
Inserting Eq. Q in Eq. pU|) , one easily finds for k = j 

for k = j ±1 



(10) 



(11) 



Q^jii) 



L-3 



dQij,k)P[Q{j,k)]Q{j,k)J2 



L-3 \ ip/ini-p/L) 



3— n 



l + Qij,k) + n 



and similarly, for k ^ j, j ± 1, 



Qij,k)\ P f ( IT T A\ 

= J ■ h[P/L,L~ A), 



h[p/L,L-i), (12) 



(13) 
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where 

n— ^ ^ 

Explicit formulas for the functions /2 and are given in the Appendix. In the limit of large L they tend to the 
simple forms 

F,{p) = hm Mp/L, L - 3) = i - ^i^f^ 



L^oc p p^ 

and 

1 2 2n — p^Pl 
F3(p) = lim h{p/L, L-4) = 2+ 3 . 

L-»oo J9 p'' 

that go to 1/2 and 1/3, respectively, in the limit p — ^ 0. 

The dynamical rule corresponding exactly to rates lO is easily found. At each time step one site (site i) of a 
one-dimensional lattice is randomly selected. Then, with probability f2{p/L, L — 3) one of the two nearest neighbors 
(i + 1 or i — 1) is chosen and a-i is set equal to it. With probability {p/ L)f^{pl L~ 4) instead, one randomly chooses 
one of the other L — 3 sites and sets <Ti equal to it. In this way, the effective topology over which the dynamics takes 
place changes at each time step. We call such network "annealed" WS network. The average properties are those of 
the original "quenched" WS network, but no permanent connection exists between non nearest neighbor sites. 

We can now write down the explicit form of the equations of motion for the mean spin and the correlation function 
on the annealed Watts- Strogatz network 

sj = -sj + f2{p/L, L-S) (s,+i + + I hip/L, L - 4) ^ (15) 

Cj.fe = —'2Cj^k + f2{p/L, L — 3) (Cj+i^fe + Cj-i^k + Cj,fe+i + Cj^k-i) + 

+ ^Mp/L,L-i)l J2 E ^^A- (16) 

\/#ij±i ;#fc,fc±i / 

For simplicity, here and in the following, the overbar is omitted. Also the arguments of /2 and /a will often be omitted. 

Eq. H16|l is complemented by the boundary condition Cj,j = 1 and by the initial condition Cj,k{t = 0). We consider 
an initial fully uncorrelated state Cj^k{t = 0) = dj^k- Hence the correlation function depends only on r — \j — k\ for 
all times and the equation of motion is 

C{r) = -2 (^1 + C(r) + 2 (^h - [C(r + 1) + C(r - 1)] + 2^ ^(0- (17) 

where the relation 2f2{p/L, L — 3) + {L — 3)p/ Lf3{p/ L, L — 4) — 1 = 0, proved in the Appendix [Eq. H5U|) ]. has been 
used. 

If we sum Eq. H15|l over j and divide by L, we obtain the temporal evolution of average total magnetization 

^ = [-1 + 2f2{p/L, L-3) + {L^3) p/L h{p/L, L - 4)] M. (18) 

Since the coefficient on the right hand side vanishes [Eq. (|50|) ] the average total magnetization is conserved, as in the 
dynamics on regular lattices. 

B. The stationary state 

We now turn to the analysis of Eq. (|17|l and consider first the continuum limit in real space, yielding 

C _2 (1 - 2A) (l + jl^) C + 2 (/, - Ijlif) C" + ^' (19) 
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Let us look for final configurations of the system, i. e. solutions of the stationary equation 

C"-AiC + $L = 0, (20) 

where 

2 _ (l~2/,)[l + 3/(£-3)] 
^""[/2-(l-2/.)/(L-3)]>' ('^^ 



and 



$L = / (22) 
^ Jo 



For finite L, taking into account that C(r = 0) = 1 and |C(r)| < 1 Vr, one obtains 



C{r) = ^ + 



exp[-ALr-]. (23) 



On the other hand, in the thermodynamic limit L —^ oo, 1/L drC{r) ~ C(oo) so that the solution is 



Imposing the consistency of Eq. H23|l with Eq. 122|l , the final correlation function turns out to be 

C(r) = 1 Vr > 0. (24) 

The final configuration of the system is a frozen fully ordered state. 

h 

C(r) = C(oo) + [1 - C(oo)] exp[-A(p)r], (25) 
with X^{p) = limi^oo A|^. From Eq. (|19() one obtains C{oo) = and since C(oo) = for t = we have 

C(r) = exp[-A(p)r]. (26) 
Hence the stationary configuration of the system is disordered, with a correlation length 

e = — = ^ = 1 (27) 

Hp) ^l/F2{p) - 2 ^py{p-l + e-P)-2' 

where the explicit form of F2{p) = lim^^oo f2ip/ L, L — 3) is computed in the Appendix. In the limit of small p the 
correlation length diverges as 

C. Preasymptotic dynamics 

We now study the equation of motion (|17|l in Fourier space by closely following the treatment for the Ising model 
on a one-dimensional lattice [19] . In this way, not only the stationary state but also the preasymptotic dynamics can 
be analyzed. Introducing 

fe' 

with k' — 27rn/L, n — ~L/2, L/2, multiplying both sides of Eq. (|17|l by e"''^'', summing over r from 1 to L — 1, we 
obtain 

^ E E = -2^Z ^ £ e'('='-'=)'- + 2i?i 5^ c,. (0(6^'=' +e'^'') e^^'^' + 2^300 E e'^'^'' 

r—1 k' r—1 k' r=l r—1 

Using Yl!^=i e'*''"''') = LSk^k' - 1, we get 

ck{t) = -7fcCfc(t) + A{t) + 2^co{t)iL6k,o - 1), (29) 
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where 



7fe =2 



1 + — - 2 I /2 - — ) cos(/c) 



and 



(30) 



(31) 



The disordered mitial condition impHes Cfe(t = 0) = 1 for all k. The boundary condition C(r = 0,t) = 1 implies 



' L 

k 



^E^'^W (32) 

k 

for all t. Eq. H29|l can be solved by Laplace transform methods. Introducing 

/•OO 

Cfc(s) = / Ck{t)e~'*dt. (33) 



Eq. (123) becomes 

scfc(s) - 1 = -7fcCfc(s) + a(s) + 2p/3Co(s)(4,o - 1/i), (34) 
where a(s) = 7fcCfc(s) is the Laplace transform of A{t). The coefficients Cfc(s) can be formally written down 

Cfc(s) = < (35) 

a(a) + l-7oco(a)/£ -z^ 

I s+7fc ' 

Notice that 70 = 2^/3. To determine a(s) we transform the boundary condition H32|) . use the expressions H35|) and 
replace the discrete sum over k with an integral (which is correct in the limit of large L we are interested in) , obtaining 

i.i<fi±i+(i-^)i„„+i,r* ' ,36, 

s is + 70 V Ls + 70 / 7-7r 27r s + jk 

For large L we can take jk{L) — > jk{L — 00) = 2[1 — (1 — 7o/2) cos(fc)]. The integral is easily carried out yielding 
[s^ +4:3 + 4[1 - (1 - 7o/2)^]]"^/2. g^j.g interested in t > 1/p ^ s <C 70 so that we obtain 



- = m + 1] 

s 



Ls + 70 2^/75' \ Ls + 70 



70 



(37) 



From Eq. H37|l . one realizes the existence of a temporal scale r ~ L/^/j^ separating two different regimes. 
For t ^ T, a{s) + 1 = 2y/jo/s so that ca{t) = and Ck{t) = 2^/jo/'^k- The correlation function is 

C{r,t) = 2^- + 2^ — ' — (38) 

L 27r 2[1 - (1 - 70/2) cos(fc)J 

The first term is negligible because t <^ t. Hence C(r) does not depend on time and decays exponentially with r. 
This is a quasi-stationary state with coexisting domains of fixed size. Such a regime lasts for a time r which diverges 
for L — > 00, so that in the thermodynamic limit the system remains asymptotically in this disordered state. 

For t ^ T instead, all Cfc(s) vanish except Co(s) = L/s. Hence C(r) = 1 and the system becomes completely ordered. 
This is the asymptotic regime for a system of finite size. 

Let us summarize the results obtained in this section. We have solved the dynamics of the voter model on an 
annealed small-world topology: in this way we have found that finite systems remain in a disordered state for a time 
proportional to their size, and then converge to the totally ordered configuration; infinite systems instead reach a 
disordered final stationary state, with a correlation function that decays exponentially over a distance ^p. For small 
values of p the correlation length diverges as 
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FIG. 1: Main; Plot of the fraction of active bonds ua for p = 0.1 and system sizes L 400, 600, 800, 1000, 1400, 2000 (from left 
to right) for the annealed case. Data are averaged over 1000 different realizations. Inset: The duration r of the plateau in the 
main part of the figure (symbols) plotted versus L. t is evaluated as the time at which ua drops below 0.01. The straight line 
is a power-law regression with best-fit exponent equal to 0.98 ± 0.03. 




r 



FIG. 2: Main: Correlation function C(r) in the quasi-stationary state for a system of size L = 10^ and p = 0.025 (symbols), 
compared with the analytical prediction, Eq. 12611 (solid line). Only one realization of the noise is considered. The agreement 
is excellent. Inset: The correlation length obtained numerically from the decay of the correlation function C(r) in the 
quasi-stationary state of annealed WS networks with L — 10^ (circles), compared with the analytical prediction, Eq. (I27II (solid 
line), and the same quantity computed for the quenched WS network (squares). 

IV. NUMERICAL RESULTS 

In order to validate the analytical results presented above we have performed numerical simulations of the voter 
model both on the annealed and on the quenched Watts-Strogatz topology. 

Figure n reports, for the annealed case, the temporal behavior of the fraction ua of active bonds, i. e. the fraction 
of nearest neighbor sites with opposite values of cr, for p = 0.1 and several values of L. After an initial decrease, 
typical of one-dimensional systems, a plateau sets in. The analytical treatment predicts such a preasymptotic regime 
to last for an interval proportional to L. The inset of Fig. ^ confirms the analytical finding. 

In the limit L ^ oo, the disordered state corresponding to the plateau becomes the asymptotic one. The analytical 
solution predicts an exponential form of the correlation function [Eq. 126() ] with the correlation length given by 
Eq. (|27|l . This analytical form is a very good approximation for the correlation function also for very large, but finite, 
values of L in the quasi-stationary regime (t < r) . This is shown numerically in the main part of Fig. [21 In the inset 
of the same figure we report the values of obtained numerically, which perfectly coincide with Eq. (|27|) . 

We have introduced the annealed version of the Watts-Strogatz network as an approximation for the WS network 
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FIG. 3; Correlation function C(r) in the quasi-stationary state for a quenched WS network of size L = 10^p = 0.025 (circles) 
and p — 0.2 (squares). Only one realization is considered. The solid lines are exponential fits. 

with quenched topology. Hence it is interesting to compare the results of the two cases to understand how well 
the approximation captures the behavior of the original system. In Ref. we have already performed a numerical 
investigation of the voter dynamics on the quenched WS topology. By comparing the results reported in Ref. and 
the theoretical approach presented here, we see that the annealed approximation correctly reproduces many of the 
important features of the original system, i. e. the existence of a regime with a disordered state before full order 
sets in and the linear dependence on L of the temporal scale r between them. Concerning the shape of the final 
correlation function in the case of infinite quenched networks, Fig. |31 shows that C{r) is exponential, another feature 
that is the same in the annealed and the quenched cases. What the annealed approximation is not able to capture is 
the quantitative dependence of the correlation length on p. This is shown in the inset of Fig. [2 where it is clear that, 
for small p, ~ in the quenched system (9j, while £,p ~ p"-^^^ in the annealed case. 



EXTENSION TO v>l 



In the previous sections we have considered an initial one-dimensional lattice with each site connected only to its 
nearest neighbors. When the number of neighbors connected to each site is > 1, one can quite easily extend the 
analytical calculations presented above. On physical grounds one expects the qualitative picture to remain the same. 
Also formulas are very similar. In general they involve, in the place of /2 and /a the functions f2u{p/ L ~ 1 — 2i>) 
and J2u+i{p / L , L — 2 — 2v). The equation for the average total magnetization becomes 



^ - [-1 + 2vf2,{plL, L-l-2v) + {L-l~2v) p/L /2.+i(p/£, L - 2 - 2z/)] M, 



(39) 



and M is again conserved due to the vanishing of the coefhcient on the right hand side [Eq. (|50|l ]. The equation for 
the correlation function C(r) becomes 



L 



L 



'-]Y.[C{r + ^) + C{r-^)] + 2P-^±^Y.C{l) 



L-l 



1=1 1=0 

Considering the continuum limit of Eq. l|4Uf) and looking for the stationary state one obtains 



where now 



and 



(1 - 2iyf2,)[l + i2iy + 1)/(L - 2i^ - 1)] 
[/2. - (1 - 2i//2.)/(i - 2j. - 1)K 

.2 i^(i^+ l)(2iy+ 1) 



>0, 
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(40) 

(41) 
(42) 

(43) 
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FIG. 4: Main: Plot of the fraction of active bonds ua for L = 3000 and p = 0.05, p = 0.025, p = 0.01 and p = 0.008 (from top 
to bottom) for the annealed case. Data are averaged over 1000 different realizations. Inset: The duration r of the plateau in 
the main part of the figure (symbols) plotted versus p. r is evaluated as the time at which ua drops below one tenth of the 
value during the plateau. The straight line is a power-law regression with best-fit exponent equal to 0.47 ± 0.03. 



Again the only solution for finite L is the completely ordered state C(r) — 1, while in the thermodynamic limit 
C(r) = exp[— A(p)r]. The correlation length is 



= 1/Xip) = 



1 



(44) 



where F2^{p) = lim^^oo f2u{p/L, L — 2v — 1). In the limit of small p, by expanding Eq. (|48|l . one finds 

~ {2vp)-'/\ (45) 

For what concerns the preasymptotic dynamics, the only formal change in Eq. (|29fl is that /s is replaced by f2u+i, 
but now the form of 7^ is different 



Ik = 2 



Pf2u+1 ^ f r Pf 2u+l 



cos(jfc) 



(46) 



By applying the Laplace transform one gets an equation formally equal to Eq. I|36|l . For generic v we cannot perform 
explicitly the integral appearing in Eq. (|36|l and hence we cannot write down the analogue of Eq. H37I) . However, 
we can guess that the only change will be the replacement of the factor l/2y^ with some other factor independent 
from L. Therefore the existence of two regimes separated by a temporal scale r proportional to L will be preserved. 
Moreover, on physical grounds, we expect the proportionality factor to scale as p~^/'^ also for generic v. To confirm 
this, we have performed numerical simulations of a system with next-nearest neighbor connections {v = 2). The 
results, presented in Fig. 01 confirm the expectation. The temporal scale r separating the quasi-stationary disordered 
state from the totally ordered configuration scales as Lp^^^^. We can conclude that the behavior of the voter model 
on the small world topology is qualitatively the same, regardless of the number v of connections between neighbors. 



VI. CONCLUSIONS 



In this paper we have studied analytically the voter model on the small-world topology. We have considered an 
annealed version of the Watts-Strogatz network, where long-range connections are not fixed, but chosen randomly at 
each time step. In this way each realization of the voter model dynamics takes place in an effective average small- world 
topology and this allows the exact solution of the equation for the correlation function of the system. 

The dynamical behavior of the model on the annealed topology is very similar to the behavior on the quenched 
WS network. Systems of finite size converge asymptotically to a totally ordered frozen state after an intermediate 
quasi-stationary stage, characterized by a finite correlation length. The duration of this preasymptotic regime is 
proportional to the number of sites. Hence systems of infinite size never reach the ordered state and remain in a 
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disordered stationary state with finite domain size. All these features are exactly the same both on the annealed and 
on the quenched versions of the network. A quantitative difference arises only in the dependence of the correlation 
length on the probability p of having a long-range connection. This discrepancy is not surprising, since a similar 
disagreement between the annealed and the quenched case has been noted previously for diffusion on Watts-Strogatz 
networks ^nd the voter model is well known to be related with first passage properties of random walkers |l8l |. 
In that case, the crossover time separating short and long time behavior of the mean number of distinct sites visited 
scales as in the quenched network and as l/p in the annealed one. This difference is due to the fact that in 
the quenched case a walker has to diffuse over a distance l/p before reaching a shortcut and deviating from the 
one-dimensional behavior. This clearly requires a time In the annealed case instead, the time needed to perform 
a long-range jump scales as l/p. In the voter model on small-world topology the boundaries between ordered domains 
perform a one-dimcnsional random walk for short times. The behavior changes when the walkers make a long range 
jump. As mentioned above, this requires a time l/p in the annealed network and p^^ in the quenched one. This 
difference generates the different scaling of the correlation length in the two types of topology. 

Despite this discrepancy the voter model dynamics on small-world networks is relatively insensitive to the quenched 
or annealed nature of the topology. An interesting question for future work is whether this insensitivity extends also 
to other ordering processes on other complex networks. 



APPENDIX 

In this Appendix we give explicit formulas for the functions /2i/(p/£, L — — 1) and /2i/+i(p/£, L — — 2). Let 
us consider 



N 



71=0 



N \ a"(l - a) 



R + n 



(47) 



For a 7^ 



fR{a,N) 



-4 r dss'^-^s + 1 - a] 
Jo 



N 



The derivation of Eq. (|48|) is easy 



R 



R + n 



n=0 



(48) 



(49) 



R+n-l 



Jo 



N 



n 



a -a) 



N — n ^71 



n=0 

dss^~\s + 1 - a)^ . 



Using Eq. l|T7jl it is easy to verify that 

Rfnia, N-R) + aiN- i?)/fl+i(a, N-R- 1)^1. 



(50) 



The vanishing of the coefficient on the right hand side of Eqs. H18|l and H39(l is obtained by setting R = 2v, a = p/L 
and iV = L - 1. 

Using Eq. H48|l the explicit formulas for /2 and /a are readily found 



/2(p/£,L-3)= ( - 



1 



(l-p/L) , {1-p/L) 



L - 1 



L-2 



L-l 



(L-2)(L-1) 



(51) 



f3{p/L,L-4)^ [ ^ 



2(1 -p/L) , {l-p/iy 



L-l 



L 



2(1 -p/L) 



L-l 



L-3 



{L - 3){L - 2){L - 1) 



(52) 
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Finally the asymptotic forms of the functions for L — > cx) are 



F2{p) = lim f2{p/L,L 



3) 



1 (1 - e^n 

p p^ 



(53) 



and 



F3(p) = lim /3(p/L,L-4) 



P 



1 



2 



+ 




(54) 



that go to 1/2 and 1/3, respectively, in the limit p — s- 0. 
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